
 
source("U:\\My Documents\\PERCENTAGE BEND R\\Rallfun-v20.R") 
library(foreign)
x <- read.spss("U:\\My Documents\\PERCENTAGE BEND R\\PERC BEND TEST DATA.SAV")

attach(x)

    x
pbcor(x$x,x$y,beta=0.2)
$cor
[1] 0.2390299

$test
[1] 1.073011

$siglevel
[1] 0.2967058

$n
[1] 21

cor(x$x,x$y)
[1] -0.03472996
jpeg(file="U:\\My Documents\\PERCENTAGE BEND R\\pbdat.jpeg")
plot(x$x,x$y)
dev.off()

